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INTRODUCTION 


The  military  Is  becoming  aware  of  Its  need  for  real-time  remote  wind 
Information  In  various  kinds  of  weather.  Wind  informaton  of  this  kind  Is 
useful  for  wind  correction  of  unguided  ballistic  projectiles  for  improved 
first  round  hit  capability. 

A  high  relative  humidity  condition- -all  atmospheric  water  essentially  in  the 
vapor  phase— is  felt  to  be  the  meteorological  condition  in  which  this  specific 
information  is  of  the  greatest  value  to  a  military  field  unit.  The  method 
chosen  to  make  the  wind  measurement  oust  be  capable  of  producing  information 
appropriate  for  the  ballistic  correction  required.  An  ungulded  ballistic 
device,  for  best  results,  requires  a  crosswind  correction  measured  in  the 
projectile  flight  path  in  real-time.  Generally,  the  first  few  hundred  meters 
is  the  critical  portion  of  the  flight  path  as  relates  to  ballistic  wind 
correction.  Methods  are  referenced  for  measuring  the  wind  velocity  remotely 
by  using  atmospheric  dust  as  a  wind  tracer. 

This  paper  presents  an  investigation,  however,  that  is  believed  to  be  the 
first  attenpt  to  use  water  vapor  as  a  remote  wind  tracer  and  discusses  the 
theory  and  feasibility  of  a  remote  wind  sensing  procedure  based  on  measuring 
the  motion  of  water  vapor  microstructui  e  by  the  Differential  Absorption  Lldar 
(DIAL)  technique.  The  ordinary  Raman  back scatter  technique  is  also 

considered.  An  example  is  given  for  a  measurement  of  the  wind  at  a  range  of 
500  m  that  appears  to  be  practical  for  a  total  measurement  time  of  under  15  s. 

BACKGROUND 

In  the  past,  wind  corrections  for  ungulded  ballistic  projectiles  were  taken  at 
a  single  point  with  propeller  anemometers.  Occasionally,  pibal  (pilot- 
balloon)  observations  of  wind  measurements  were  used  for  wind  correction  above 
the  surface.  For  long  range  projectiles,  radiosonde  wind  information  was  used 
at  higher  altitudes.  These  various  measurements  lagged  badly  in  time,  and  the 
measurement  position  was  poor.  The  fact  that  these  measurement  types  are  not 
taken  in  the  path  that  the  projectile  will  travel  Is  still  another 
disadvantage.  A  measuring  technique  that  will  be  far  more  appropriate  for  the 
best  correction  Is  one  that  will  give  real-time  Information  of  the  wind 
profile  along  the  projectile's  path.  A  remote  wind  measuring  method  can 
fulfill  these  criteria  much  better  than  past  techniques. 

In  the  past  several  years,  researchers  have  reported  several  single  ended 
remote  Ildar  wind  measuring  techniques  In  the  open  literature.  Schwlesow  and 
Cupp1  and  Post  et  al 1  demonstrated  the  feasibility  of  using  the  Doppler 
technique,  with  a  continuous  wave  (cw)  CO*  lldar,  for  measuring  the  horizontal 


‘Schwlesow,  R.  L.,  and  R.  E.  Cupp,  "Remote  Doppler  Velocity  Measurements  of 
Atmospheric  Dust  Devil  Varieties,"  Appl  Opt,  15:1,  1976. 

*Post,  M.  J.,  et  al,  "A  Comparison  of  Anemometer  and  Lldar  Sensed  Wind 
Velocity  Data,"  J  Appl  Meteorol,  17:1179-1181,  1978. 
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wind  component  parallel  to  the  laser  beam.  With  this  specific  method, 
however,  the  most  Important  component  in  ballistic  wind  correction— the  wind 
component  perpendicular  to  the  laser  beam— cannot  be  measured. 

Teoste  and  Capes,*  using  the  Doppler  technique  with  a  pulsed  CO*  lidar  that 
conically  scans,  showed  the  potential  of  such  a  device  for  producing  a 
vertical  wind  profile  of  the  horizontal  components  to  an  altitude  of  8  to  10 
km.  This  approach  is  limited  because  it  will  obtain  only  the  wind  profile 
directly  over  the  device,  not  the  perpendicular  wind  component  across  the  path 
of  interest.  Schotland*  developed  a  Doppler  visible  light  Ildar  in  which  an 
argon  laser  with  a  modulated  cw  beam  is  used  to  observe  Doppler  shifts  in  the 
modulation  frequency.  This  approach,  too,  is  limited  because  it  will  measure 
only  the  velocity  vector  along  the  laser  beam. 

Kunkel  et  al  *  and  Sroga  et  al  *  reported  on  a  technique  for  measuring  wind 
velocity  and  direction  in  a  horizontal  plane  by  tracking  tenuous  natural  dust 
clouds  in  the  atmosphere.  Using  a  scanning  lidar,  this  approach  has  been 
demonstrated  out  to  a  range  of  8  km. 

Derr  et  al 7  discussed  the  feasibility  of  the  transit-time  lidar  (TTL) 
technique  and  demonstrated  the  principle  with  some  pronounced  dust 
irregularities.  But  these  techniques  appeared  to  be  slow  because  they  depend 
upon  observing  the  scattering  patterns  from  the  dust  in  two  specific 
volumes.  Armstrong  et  al •  and  Barber  et  al*  examined  the  TTL  method  as  a  wind 
measuring  technique  for  accuracy  in  an  enhanced  dust  condition.  Performed  at 
a  range  of  250  m,  this  research  discussed  sample  rates  and  volume  spacings  as 
well.  The  TTL  method  does  produce  the  wind  vector  perpendicular  to  the  laser 


*Teoste,  R.,  and  R.  N.  Capes,  "High-Altitude  Infrared  Radar  Wind 

Measurements,"  J  Appl  Meteorol,  17:1575,  1978. 

"Schotland,  R.  M.,  "The  Duel  Frequency  Doppler  Lidar  Technique  for  Wind 
Measurement,"  paper  presented  at  the  9th  International  Laser  Radar  Conference, 
Munich,  Germany,  1979. 

*Kunkel,  K.  E. ,  E.  W.  Eloranta,  and  J.  A.  Weinman,  "Remote  Determination  of 
Winds,  Turbulence  Spectra,  and  Energy  Dissipation  Rates  in  the  Boundary  Layer 
from  Lidar  Measurements,"  J  Atmos  Sci,  37:978-985,  1980. 

•Sroga,  J.  T.,  E.  W.  Eloranta,  and  T.  L.  Barber,  "Lidar  Measurement  of  Wind 
Velocity  Profiles  In  the  Boundary  Layer,"  J  Appl  Meteorol ,  Vol  19,  No  5,  1980. 

7Derr,  V.  E.,  and  C.  G.  Little,  "A  Conparison  of  Remote  Sensing  of  the  Clear 
Atmosphere  by  Optical,  Radio,  and  Acoustic  Radar  Techniques,"  Appl  Opt, 
9:1976,  1970. 

•Armstrong,  R.  L.,  J.  B.  Mason,  and  T.  L.  Barber,  "Detection  of  Atmospheric 
Aerosol  Flow  Using  a  Transit-Time  Lidar  Velocimeter,"  Appl  Opt,  15:2891,  1976. 

’Barber,  T.  L.,  and  J.  B.  Mason,  "A  Transit-Time  Lidar  Wind  Measurement,  A 
Feasibility  Study,"  EC0M-5550,  Atmospheric  Sciences  Laboratory,  US  Army 
Electronics  Command,  White  Sands  Missile  Range,  NM,  1974. 


beam.  Barber  and  Rodriquez1*  demonstrated  the  capability  of  the  TTL  In 
extremely  clean  air  conditions.  With  this  technique,  the  measurement  point 
can  be  placed  at  the  point  of  interest,  and  several  measurements  along  the 
laser  beam  can  be  taken  simultaneously.  It  is  quite  applicable  for  ballistic 
projectile  wind  correction,  because  the  sampling  period  is  relatively  short 
(10  to  20  s),  and  the  sampling  point  can  be  positioned  in  the  path  of 
interest.  All  of  these  remote  sounding  methods  depend  on  tracking  the  flow  of 
natural  dust  or  haze  to  measure  the  wind.  This  report,  however,  considers  the 
practicality  of  using  water  vapor  in  a  wet  atmosphere  as  a  wind  tracer.  Some 
of  the  necessary  data  handling  procedures  for  this  method  and  for  those  using 
dust  are  similar. 


PROCEDURE 

This  paper  considers  methods  for  using  variation  in  absolute  water  vapor 
content11  12  as  a  tracer  to  remotely  measure  wind  speed  and  wind  direction  in 
a  wet  atmosphere  (relative  humidity  -  (RH)  high)- -where  water  is  essentially 
in  the  gas  phase.  Investigators  have  remotely  measured  water  vapor  content  11 
12  and  noted  microstructure  in  the  concentration.  When  a  relatively  reliable 
measurement  of  the  water  vapor  content  can  be  made  remotely,  the 
microstructure  in  the  water  vapor  can  be  obtained. 

Measurement  of  the  microstructure  will  be  done  at  two  well-defined  fixed 
positions  in  the  atmosphere  a  distance  "D"  apart  (figure  1).  Various 
mathematical  techniques  can  be  used  to  obtain  the  time  for  the  microstructure 
to  traverse  "D."  With  the  distance  and  time,  the  wind  speed  and  wind 
direction  perpendicular  to  the  laser  beam  can  be  derived. 

Two  possible  techniques  for  sensing  the  water  vapor  content  at  a  point  in  the 
atmosphere  were  considered:  (1)  ordinary  Raman  scattering  technique11  and  (2) 
DIAL.12  “ 

Although  the  ordinary  Raman  scattering  process  will  produce  the  necessary 
information,  the  inefficiency  of  this  process  requires  considerable  laser 
power.  This  necessary  power  level  creates  problems  for  eye  safety  and 
equipment  size.  Investigators  using  the  DIAL  technique  have  found  that 


‘•Barber,  T.  L. ,  and  R.  Rodriquez,  "Transit-Time  Lidar  Measurement  of  Near 
Surface  Winds  in  the  Atmosphere,"  ASL-TR-0033,  US  Amy  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  1979. 

“Cooney,  John  A.,  "Comparisons  of  Water  Vapor  Profiles  Obtained  by  Radiosonde 
and  Laser  Backscatter,"  J  Appl  Meteorol ,  10:301,  1971. 

“Schotland,  R.  M.,  "Errors  in  the  Lidar  Measurement  of  Atmospheric  Gases  by 
Differential  Absorption,"  J  Appl  Meteorol ,  13:71,  1974. 

“Wright,  M.  L.,  et  al,  "A  Preliminary  Study  of  Air  Pollution  Measurement  by 
Active  Remote  Sensing  Techniques,"  Final  Report,  Stanford  Research  Institute, 
Menlo  Park,  CA,  1966. 
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Figure  1.  Sample  volumes.  A  graphic  portrayal  of  the  scattering  volumes  in 
the  atmosphere.  Both  wavelengths  A  and  B  will  be  transmitted  along 
path  E.  The  lidar  pointing  angle  will  be  changed  by  some  small 
amount  &e,  moving  to  E'.  Wavelengths  A  and  B  will  be  transmitted 
again  along  E' .  The  position  along  the  laser  beam  can  be  adjusted 
by  the  time  delay  In  the  gate  circuit.  0,  the  beam  spacing,  can  be 
calculated  by  D  s  R  sin  40.  The  total  system  is  shown  in  figure  2. 
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measuring  water  vapor  in  the  atmosphere  appears  to  be  the  more  promising 
method.12  12  l*  14  14  This  technique  involves  transmitting  two  pulsed 
wavelengths  of  light  along  a  common  path:  the  first  wavelength  of  light,  (A), 
is  essentially  not  absorbed  by  water  vapor;  the  other,  (B),  is  appreciably 
absorbed  by  water  vapor. 

Preferably,  the  wavelength  B  should  be  exactly  centered  on  a  water  vapor 
absorption  line.  The  lidar  receives  backscatter  from  the  dust  and  haze  at 
both  wavelengths  A  and  B  (figure  1  and  figure  2,  E)  from  two  common  scattering 
volumes,  F  and  6,  making  it  possible  to  measure  the  water  vapor  content  in  the 
volume  of  interest,  aR  in  the  interval  from  F  to  G.  For  best  results, 
wavelengths  A  and  B  should  not  be  further  apart  in  wavelength  than  about  0.3 
cm"1,12  so  that  the  volume  backscatter  coefficient  from  the  dust  and  haze  in 
the  atmosphere  in  volumes  F  and  G  will  essentially  be  equal.  The  laser  that 
transmits  wavelengths  A  and  B  must  have  a  pulse  length  capable  of  sufficient 
range  resolution  for  the  approach  discussed  here,  less  than  30  ns  pulse 
length.  Both  wavelengths  A  and  B  must  be  transmitted  along  the  same  path  in 
the  atmosphere  (figure  1,  E).  These  pulses  must  be  no  more  than  a  few  hundred 
microseconds  from  each  other*  to  minimize  the  problem  of  changes  in 
atmospheric  particulate  concentration  with  time  in  the  scattering  volumes. 
The  lidar  receives  back  scattered  light  from  scattering  volumes  F  and  G  at  both 
wavelengths  A  and  B.  These  returns  will  be  called  PAfr  and  PAG  for  wavelength 

A,  and  Pqjt  and  P0G  for  wavelength  B.  When  transmitting  wavelength  A,  the 

signal  received  from  the  2  volumes,  F  and  G,  will  be  essentially  scatter  from 
the  dust  and  haze  affected  only  by  the  1/R 2  term  and  scattering  losses.  R  is 
the  range  from  lidar  to  scattering  volume  F.  The  signal  received  from  the  2 
scattering  volumes,  F  and  G,  will  be  absorbed  by  water  vapor  at  wavelength 

B.  Considering  one  specific  beam  path  (figure  1,  E)  with  a  pulse  A  and  a 
pulse  B,  four  specified  scattered  returns  will  be  obtained--PAF,  PAG,  PBF,  and 


12Schotland,  R.  M.,  "Errors  in  the  Lidar  Measurement  of  Atmospheric  Gases  by 
Differential  Absorption,"  J  Appl  Meteorol,  13:71,  1974. 

14Wright,  M.  L. ,  et  al ,  "A  Preliminary  Study  of  Air  Pollution  Measurements  by 
Active  Remote  Sensing  Techniques,"  Final  Report,  Stanford  Research  Institute, 
Menlo  Park,  CA,  1966. 

1%Fraley,  P.  E. ,  and  K.  N.  Arallal  Rao,  "High  Resolution  Infrared  Spectra  of 
Water  Vapor  V,  and  V*  Bands  of  H^'O,"  J  Mol  Spectroscopy,  29:312-316,  1969. 

**Hinkley,  E.  D. ,  ed.  Laser  Monitoring  of  the  Atmosphere,  (Berlin,  Heidelberg, 
New  York:  Springer-Verlag,  1976). 

‘•Browell,  E.  V.,  T.  D.  Wilkerson,  and  T.  J.  Mcllrath,  "Water  Vapor 
Differential  Absorption  Lidar  Development  and  Evaluation,"  Appl  Opt,  18:3474, 
1979. 

’Barber  T.  L. ,  and  J.  B.  Mason,  "A  Transit-Time  Lidar  Wind  Measurement,  A 
Feasibility  Study,"  EC0M-5550,  Atmospheric  Sciences  Laboratory,  US  Army 
Electronics  Command,  White  Sands  Missile  Range,  NM,  1974. 


9 


PBq.  By  using  the  DIAL  technique,  a  value  for  the  abso'.ute  water  vapor 
content  in  &R  can  be  derived. 

THEORY 

To  disucss  the  fundamentals  of  remotely  measuring  the  water  vapor  content, 
using  the  DIAL  technique,  the  first  step  is  to  consider  the  fundamental  lidar 
equation: 


P  =  PQ  s(R)XR-2  exp  -2  /  a(R)dR 


Where  P  is  the  received  power  by  the  lidar,  PQ  is  the  transmitted  power,  C  is 

the  speed  of  light,  and  T  is  the  laser  pulse  length  in  time.  b  is  the 
backscatter  coefficient  of  the  atmosphere,  X  is  the  effective  frontal  area  of 
the  lidar  receiver,  and  R  is  the  range  to  the  scattering  volume  of  interest. 

D 

The  term  -2  /  a(R)dR  gives  the  extinction  along  the  path  R  in  the 

atmosphere.  0  a  is  the  unit  volume  extinction  coefficient  in  a  specific 
volume  of  the  atmosphere. 

The  second  step  is  to  break  down  the  unit  volume  extinction  coefficient  at 
wavelength  B  into  two  parts: 


°B  =  aS  +  °w 
and 

°A  =  °S  . 


(2) 


Wavelengths  A  and  B  have  been  previously  defined.  aw  is  the  unit  volume 

extinction  coefficient  at  B  caused  by  water  vapor  absorption.  The  a$  is  the 

unit  volume  extinction  coefficient  at  wavelengths  A  and  B  caused  by 
particulate  and  molecular  scattering  and  minor  absorption  by  materials  other 
than  water  vapor.  This  statement  concerning  a$  is  appropriate  since  A  and  B 

are  no  further  apart  in  wavelength  than  0.3  cm"1. 


The  third  step  is  to  find  the  relationship  between  total  water  vapor  and  aw: 


w 


Ww 


(3) 


n 


Kw  is  the  water  vapor  absorption  cross  section  at  wavelength  B.  ^  is  the 
total  water  vapor  content  in  1  unit  volume. 

Next  set: 


and 


CT 
AO  T 

CT 
BO  7~ 


X 

X  . 


(4) 


PA0  is  output  power  of  the  lidar  at  wavelength  A,  and  PB0  is  the  output  power 
of  the  lidar  at  B. 

With  the  criteria  for  A  and  B  as  previously  stated, 


Ba(R)  =  BB(  R)  =  b(R)  . 


(5) 


Next,  two  equations  are  written  giving  the  return  at  A  and  B, 


PAF  =  CA 


6(R)  exn 
~JT  eXp 


m  m 

-2  ;R  a<«  (R)dR  . 
0  5 


/R  CaS(R)  +  °w(R,]dR|* 


(6) 

(7) 


The  terms  Z(R)  and  Z(R  +  DR)  are  now  defined  as. 


and 


Z(R)  *  In  PBF  -  In  PAF  and  Z(R  +  aR)  =  In  Pgg  -  In  PAG 
Z(R  +  DR)  =  In  Pgg  -  In  PAQ  . 


(8) 


Inserting  (6)  and  (7)  into  (8)  the  resulting  two  equations  are: 
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1 


Z(R) 


and 


Z(R  +  aR)  =  1 


'  ’"W' 

n(^)' 


2  IK  cr  (R)dR  , 
o 


2  /(R+AR)a  (R)dR 
0  w 


(9) 


If  measurements  are  made  at  two  points  (figure  1,  F  and  G),  at  a  range  R  and 
range  (R  +•  aR),  to  obtain  the  difference  over  the  length  aR  becomes: 


Z(R  +  aR'  -  Z(R)  =  -2  /R+aRo  (R)dR  . 

n  W 


(10) 


The  NwaR  (equation  (3))  can  be  obtained  by: 


.R+aR 


v"  -  K«  *  C  vR'dR  • 

K 


(ii) 


The  aw  is  the  average  value  for  ow  over  aR.  Equation  (11)  becomes: 


N»SR  ■  r;  '««  • 


(12) 


nwaR  ‘is  the  total  water  vapor  in  the  interval  aR.  Rewriting  equation  10 
results  in: 


or 


Z(R)  -  Z(R  +  aR)  =  2  ;(R+aR)0  ( R)dR 

R  w 


—  ~  *(R  *  =  /(R+AR£  (R)dR 


<jwaR 
NwaR  Kw 


(13) 


(14) 
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Note  that  the  right  portions  of  equations  (11)  and  (14)  are  equal.  Multiplied 
by  2/Kw  equation  (14)  becomes: 


Z(R)  -  Z ( R  +  a)  _  „„ 

IT  *NwaR  (15) 


Insert  the  values  for  Z(R)  and  Z(R  +  aR)  from  equation  (8),  then 


l  =  1 

'war 


(16) 


or 


I  -  1 

'war  2kw 


(17) 


Rearranged,  this  equation  becomes: 


NwaR  =  FT 
w 


AR, 


and 


"wsR  ■  °w  lR 


Previously  defined,  these  four  values  for  P  are  the  four  signal  values 
measured  with  a  lidar.  As  formerly  stated,  the  N^r  is  the  absolute  quantity 

of  the  water  vapor  in  the  interval  aR.  This  derivation  produces  a  single 
value  of  N^r  for  each  single  set  of  PAp,  PAG,  PBp,  and  PBS. 

CALCULATIONS 

The  previous  section  gives  an  outline  of  the  theory  related  to  the  DIAL 
technique.  This  section  will  give  an  example  of  a  wavelength  region  that  can 
be  used  for  lidar  measurements  of  water  vapor.  Water  vapor  has  many 
absorption  bands,1*  and  the  band  between  715  and  732  nm  is  good  for  water 


‘♦Fraley,  P.  E.,  and  K.  N.  Arallal  Rao,  "High  Resolution  Infrared  Spectra  of 
Water  Vapor  V,  and  V*  Bands  of  H1I8O,"  J  Mol  Spectroscopy,  29:312-316,  1969. 
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vapor  measurements.1*  17  There  are  62  absorption  lines  in  this  area, 
carefully  measured  and  described  by  Wilkerson,17  that  are  appropriate  to  use 
in  the  DIAL  technique  for  measuring  water  vapor.  Because  the  individual  line 
strengths  vary  through  a  factor  of  about  150,  different  lines  can  be  used  in 
different  meteorological  conditions.  There  are  also  small  windows  in  this 
region  on  which  to  set  wavelength  A. 

The  specific  example  that  will  be  discussed  for  wavelength  B  is  the  line  at 
723.61  nm.  The  temperature  chosen  is  20°C  and  an  RH  of  90  percent,  the 
partial  pressure  at  this  temperature  and  RH  is  15.79  mm  for  the  water  vapor. 
From  Wilkerson  et  al,17  the  value  given  for  the  cross  section  of  water  vapor 
per  molecule  is  48.2  10” *  *  cm*.  With  this  value  and  the  conditions  previously 
stated,  a  unit  volume  extinction  due  solely  to  water  vapor  absorption  can  be 
calculated.  The  number  of  water  molecules  in  a  cubic  meter  of  the  atmosphere 
at  this  temperature  and  1  atmosphere  pressure  is: 


no. 


15,79  mm  273K 

760  mm  293K 


i  » 

10  /m 

- j - 

22400  cm  /mole 


.6.023  x  10  =  5.21  x  10*  molecules. 

(18) 


The  cross  section  per  molecule  at  line  center  obtained  from  Wilkerson17  and 
multiplied  by  the  number  of  molecules  in  the  specified  volume  of  interest 
gives  the  volume  cross  section: 


2  9  * 

o  *  5.21  x  10  molecules  — ,  x  48.2  10"*'cm2  =  2.51  x  lO'Vm 
w  10  cm 


(19) 


The  extinction  value  2.51  x  l0“,^^,,  falls  well  within  the  values  that  appear  in 
table  1,  column  1,  row  5. 

Another  example  that  has  been  worked  out  in  detail  in  the  literature1*  for  CO 
as  the  molecule  and  aR  of  15  m  points  to  the  relative  accuracy  of  this 
approach. 


“Browell,  E.  V.,  T.  D.  Wilkerson,  and  T.  J.  Mcllrath,  “Water  Vapor 
Differential  Absorption  Lldar  Development  and  Evaluation,"  Appl  Opt,  18:3474, 
1979. 

17Wilkerson,  T.  D.,  G.  Schwemmer,  and  G.  Gentry,  “Intensities  by  Ni  Collision- 
Broadening  Coefficients  Measured  for  Selected  Ht0  Absorption  Lines  Between  715 
and  732  nm,"  J  Quant  Spectroscopy  Radi  at  Transfer,  22:315-331,  1979. 

MByer,  R.  L. ,  "Review  Remote  Air  Pollution  Measurement,"  Optical  and  Quantum 
Electronics,  7:147-177,  1975. 
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This  case  gives  an  example  of  the  pulse  power  necessary  to  achieve  the  degree 
of  accuracy  required  in  the  region  of  10  to  100  m3  per  pulse.1'  As  an 
example,  the  range  R  will  be  chosen  as  490  m  and  aR  10  m.  The  absolute 
humidity  is  the  quantity  to  be  measured,  so  the  temperature  will  be  held 
constant  at  20°C.  An  estimate  of  the  relative  data  accuracy  is  shown  in  this 
example  from  the  portion  Inside  the  large  brackets  in  equation  (17).*'  "Given 
typical  detection  and  digitization  equipment,  a  reasonable  value  for  Ain  is 
about  0.02."15 

Several  investigators  have  examined  the  types  of  error  that  can  arise  in  such 
a  measurement. 11  ,s  *•  **  *'  In  the  same  general  type  of  measurement  as 
discussed  here,  an  absolute  accuracy  of  1  to  6  percent  is  shown  to  be 
possible.1'  "  The  portion  of  equation  (17)  inside  the  large  square  bracket 
can  be  simplified  by  considering  equation  (1).  In  equation  (17),  the  portion 
outside  the  large  square  brackets  (l/2kw)  is  a  constant  for  a  specific 

wavelength  of  B;  hence,  for  a  series  of  measurements  where  B  is  not  changed, 
the  relative  accuracy  is  not  affected.  The  first  term  inside  the  large 
brackets  in  equation  (17),  In  PAG/PA p,  can  be  restated  as: 


In 


(R 


+  ART 
^ - 


(20) 


The  second  portion  of  equation  (17)  Inside  large  brackets  can  be  rewritten  as 


i.  (R  +  AR)'2  2VR 

I  n  "  n  —  •  I  t 

R 


(21) 


or  combining  the  terms  inside  the  large  brackets  in  equation  (17): 


» 'Byer  R.  L. ,  and  Max  Garbuny,  "Pollutant  Detection  by  Absorption  Using  Mie 
Scattering  and  Topographic  Targets  as  Retroreflectors,"  Appl  Opt,  12:901, 
1973. 

1  'Hi nkl ey  E.  D.,  ed.  Laser  Monitoring  of  the  Atmosphere,  (Berlin,  Heidelberg, 
New  York:  Springer  Verlag,  1976). 

1  'Schotland,  R.  M.,  "Errors  In  the  Ildar  Measurement  of  Atmospheric  Gases  by 
Differential  Absorption,"  J  Appl  Meteorol,  13:71,  1974. 

•'Byer,  R.  L. ,  "Review  Remote  Air  Pollution  Measurement,"  Optical  and  Quantum 
Electronics,  7:147-177,  1975. 

"Ahmed,  Samir  A.,  "Molecular  Air  Pollution  Monitoring  by  Dye  Laser 
Measurement  of  Differential  Absorption  of  Atmospheric  Elastic  Backscatter," 
Appl  Opt,  12:901,  1973. 
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In 


AG 

AF 


In 


(R  VR)  !  2o  aR  =  Ain  =  In  2a  aR 

"  W 

(22) 


From  a  simple  examination  of  equation  (17),  one  can  see  that  the  total  water 
vapor  content  In  the  Interval  aR  Is  a  constant  times  Ain.  -for  a  fixed  aR, 

equation  (22)  shows  the  dependence  of  in  the  sample  volume  on  Ain. 

Equation  (17)  shows  that  for  a  given  aR  the  total  water  vapor  content  In  the 

interval  aR  is  a  constant  times  a  .  This  is  the  H  measured  in  aR.  For  best 

W  W 

results,  keeping  a  reasonable  data  quality  in  this  measurement,  wavelength  B 

should  be  chosen  so  that  ow  in  aR  will  be  maintained  in  a  preferred  range.  In 

the  example  considered  in  table  1,  aR  =  10  m  and  R  =  490  m.  As  previously 
stated,  to  maintain  a  reasonable  data  quality,  physically  reasonable  value  for 
Ain  is  0.02  or  greater.  Recall  that  different  values  of  aw  can  be  attained 

with  a  constant  water  vapor  content  in  the  sample  volume,  by  tuning  wavelength 
B  to  various  water  vapor  absorption  lines.  Six  example  water  vapor  absorption 
lines  are  given  in  table  1.  Information  in  columns  1  and  2  were  obtained  from 
Wilkerson.17  In  row  1  of  the  table,  columns  4,  5  and  6  are  significantly 
below  0.02,  making  this  wavelength  choice  an  unfavorable  one.  At  the  other 
extreme,  the  value  In  column  3,  row  6,  presents  a  different  problem.  This 
number  is  the  transmission  of  wavelength  B  over  a  range  2R.  For  the  case  as 
stated  in  column  4,  row  6,  the  measurable  signal  at  wavelength  B  will  be 
reduced  by  water  vapor  absorption  by  a  factor  of  about  23.  This  reduces  the 
measurement  reliability  of  B.  For  optimum  results,  using  these  two  effects  as 
outer  bounds,  the  preferred  wavelength  chosen  should  be  in  the  region  of  rows 
2,  3,  or  4. 


TABLE  1.  A  COMPARISON  OF  WATER  VAPOR  EXTINCTION  AND  ALn 


Wavelength 

B(NM) 

aw 

M"  1  *.l 

Tn 

Ain 

RH  65% 

Ain 

RH  90% 

Ain 

RH  115% 

1. 

717.59 

0.0045 

0.641 

0.0064 

0.009 

0.011 

2. 

717.76 

0.0109 

0.336 

0.015 

0.021 

0.027 

3. 

724.08 

0.0158 

0.206 

0.022 

0.031 

0.039 

4. 

719.39 

0.0213 

0.119 

0.031 

0.043 

0.055 

5. 

723.61 

0.0251 

0.0812 

0.036 

0.05 

0.064 

6. 

718.63 

0.0316 

0.042 

0.045 

0.063 

0.081 

l7W11kerson,  T.  D.,  G.  Schwemmer,  and  G.  Gentry,  "Intensities  by  N*  Collision- 
Broadening  Coefficients  Measured  for  Selected  Ha0  Absorption  Lines  Between  715 
and  732  nm,"  J  Quant  Spectroscopy  Radi  at  Transfer,  22:315-331,  1979. 
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Wavelength  is  the  line  center  of  a  specific  absorption  line  for  water  vapor. 
ow  is  the  extinction  in  a  1  hi  path  due  to  water  vapor  absorption  at  wavelength 
B.  TQ  is  the  optical  transmission  over  a  path  2R  at  wavelength  B,  considering 
only  water  vapor  absorption.  Ain  is  described  in  equation  (22). 

At  this  point  the  question  whether  sufficient  water  vapor  variation  naturally 
occurs  in  the  atmosphere  to  make  this  measurement  possible  should  be 

considered.  Several  investigators  have  measured  the  high  frequency  variation 
of  the  humidity  in  the  natural  atmosphere.  Some  data  have  been  taken  over 
land  surfaces 11  and  rather  extensive  data  at  different  heights  have  been  taken 
over  the  ocean.1*  11  There  is  10  to  15  percent  relative  variation,  conpared 
to  the  average  cell  size  ranges  from  5  to  50  m.  This  10  to  15  percent  should 
be  adequate  to  make  this  wind  measuring  technique  possible. 

DATA  ANALYSIS 

The  magnitude  of  the  variations  in  the  water  vapor  microstructure  mist  be  some 
minimum  value  or  greater  to  make  it  practical  to  measure  the  wind  in  this 
manner.  Remember  that  a  portrayal  of  the  water  vapor  microstructure  is 

produced  when  a  series  of  single  point  -neasurements  are  taken  in  a  moving 

medium.  With  these  various  results  in  mind,  10  percent  is  a  reasonable  value 

to  require  as  a  minimum  in  the  microstructure  variation  for  a  practical 
measurement.  Individuals  have  observed  microstructure  in  the  water  vapor  by 
remote  sounding.11  11  As  stated  in  the  last  section,  fixed  point  measurements 
have  measured  variations  approaching  20  percent  in  the  water  vapor. 

Schotland11  has  made  error  analyses  in  several  areas  of  this  technique.  The 
change  in  atmospheric  dust  with  time  Is  a  problem.  For  a  single  measurement, 
pulses  at  wavelength  A  and  B  can  be  no  further  apart  in  time  than  1  ys  to  hold 
the  error  to  10  percent.  Uncertainties  in  pressure  and  temperature  in  the 
sampling  volume  cause  an  error  of  approximately  1  percent  in  the  calculated 
absorption  cross  section.  Schotland  performed  an  error  analysis,  using  a 
tuneable  ruby  lidar  operating  at  694.38  nm,  by  comparing  the  water  vapor  line 


1  Champagne,  F.  H.,  C.  A.  Friehe,  and  J.  C.  LaRue,  "Flux  Measurements,  Flux 
Estimation  Techniques,  and  Fine-Scale  Turbulence  Measurements  in  the  Unstable 
Surface  Layer  Over  Land,"  J  Atmos  Sci,  Yol  34,  1977. 

“Phelps,  G.  T.,  and  S.  Pond,  "Spectra  of  the  Temperature  and  Humidity 
Fluctuations  and  of  the  Fluxes  of  Moisture  and  Sensible  Heat  in  the  Marine 
Boundary  Layer,"  J  Atmos  Sci,  Vol  28,  1977. 

“Pond,  S.,  et  al,  "Measurements  of  the  Turbulent  Fluxes  of  Momentum,  Moisture 
and  Sensible  Heat  Over  the  Ocean,"  J  Atmos  Sci,  Vol  34,  1977. 

“Cooney,  John  A.,  "Comparisons  of  Water  Yapor  Profiles  Obtained  by  Radiosonde 
and  Laser  Backscatter,"  J  Appl  Meteorol,  10:301,  1971. 

“Schotland,  R.  M.,  "Errors  in  the  Lidar  Measurement  of  Atmospheric  Gases  by 
Differential  Absorption,"  J  Appl  Meteorol,  13:71,  1974. 
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center  and  laser  output  for  coincidence.  The  error  estimate  varied  from  3  to 
9  percent  according  to  the  altitude  above  sea  level  where  measurements  were 
taken  for  aR  =  200  m.  With  these  various  results  in  mind,  a  reasonable  value 
to  require  as  a  minimum  in  the  microstructure  variation  is  20  percent  for  a 
practical  measurement.  Schotland  (personal  communication,  1981),  using  a  DIAL 
system  to  measure  water  vapor,  estimated  a  microstructure  variation  of  at 
least  10  percent  over  noise.  Cooney  (personal  communication,  1981),  using  a 
Raman  lidar,  had  observed  microstructure  variations  of  about  25  percent  in  the 
water  vapor  content.  In  cases  such  as  these,  it  is  possible  to  use  the  water 

vapor  variation  as  a  tracer  to  remotely  measure  wind  velocity. 

To  use  variations  in  water  vapor  content  in  the  atmosphere  to  measure  a 
component  of  the  wind  velocity  remotely  perpendicular  to  the  lidar  beam,  the 
water  vapor  variation  must  be  tracked.  A  method  for  doing  this  tracking  in  a 
moving  medium— the  wind  moving  the  atmosphere— can  be  done  by  measuring  the 
microstructure  at  two  fixed  positions  a  known  distance  apart.  The  time  of 

flight,  then,  can  be  determined  across  these  two  positions.  This  is  done  by 

first  taking  a  set  of  single  measurements  at  a  known  rate  at  2  fixed  points  in 
space.  By  plotting  the  relative  magnitude  of  these  values,  a  portrayal  of  the 
water  vapor  microstructure  can  be  obtained. 

The  time  of  flight  is  obtained  by  cross-correlating  these  two  sets  of  data, 
hence  deriving  a  time  aT,  over  the  distance  D.  The  wind  velocity,  on  a  line 
through  the  center  of  these  two  points,  is  (D/aT)  =  VD.  Several  parameters 
must  be  considered  carefully  before  taking  a  meaningful  data  set,  so  that  the 
necessary  information  to  derive  the  wind  velocity  is  available.  The  first 
question  asked  is  how  many  individual  samples  per  second  should  be  taken  for 
good  results.  This  has  been  discussed  previously*  “  for  tracking  natural 
atmospheric  dust.  A  sampling  rate  adequate  for  this  purpose  is  1000  samples 
per  second.  When  these  values  are  collected,  the  group  or  data  block  from  one 
scattering  volume  should  be  broken  down  into  smaller  groups:  for  example,  25 
individual  values.  These  values  should  be  averaged  with  each  of  the  average 
values  grouped  as  a  second  data  block  reducing  occasional  sharp  irregularities 
in  such  data.  These  irregularities11'  are  felt  to  be  caused  by  optical 
turbulence  in  the  laser  beam,  from  lidar  to  sample  volume.  With  these  2  new 
blocks  of  average  data  points,  one  from  each  of  the  2  fixed  positions,  the 
next  step  to  consider  is  how  much  total  time  the  block  should  cover.  Here  it 
might  be  imagined  that  the  longer  the  total  period,  the  easier  to  determine 
aT.  This  is  not  the  case  because  the  wind  is  a  constantly  changing 
parameter.  Considering  this,  along  with  past  experience  in  similar  lidar 
techniques,  a  total  time  of  about  10  s  is  reasonable.  If  during  this  period 
of  time,  however,  the  wind  velocity  changes  radically  when  the  two  data  blocks 
are  cross-correlated,  this  change  substantially  degrades  the  cross-correlation 
results. 


’Barber  T.  L.,  and  J.  B.  Mason,  "A  Transit-Time  Lidar  Wind  Measurement,  A 
Feasibility  Stu<fy,"  EC0M-5550,  Atmospheric  Sciences  Laboratory,  US  Army 
Electronics  and  Development  Command,  White  Sands  Missile  Range,  NM,  1979. 

“Barber  T.  L.,  and  R.  Rodriquez,  "Transit-Time  Lidar  Measurement  of  Near 
Surface  Winds  in  the  Atmosphere,"  ASL-TR-0033,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  1979. 
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Another  consideration  in  this  methodology  is  the  value  for  D  (figure  1). 
There  are  two  limits.  If  D  is  too  small,  for  example  1  m,  this  limits  the 
higher  wind  velocities  that  can  be  measured.  With  the  previously  stated 
average  value  rate  25/1000  =  1/40  s,  the  two  highest  wind  values  that  are 
measurable  are  40  m/s  and  20  m/s. 

If  a  value  for  D  is  too  large,  for  example  300  m,  the  microstructure  can 
change  appreciably  from  one  volume  to  the  next.  "0"  should  occur  in  the  range 
5  to  15  m,  generally  the  same  as  aR.  The  total  number  of  points  that  most 
computer  systems  will  accept  In  the  cross-correlation  procedure  adds  further 
difficulty  in  this  area.  This  procedure  generally  is  4096  points  In  each  of 
the  two  averaged  blocks.  For  each  averaged  point  in  the  example  given  above, 
1/40  s  period  is  covered  or  40/s.  In  10  s  of  time,  this  will  be  400  values. 

In  most  cross-correlation  procedures,  512  points  can  be  used.  Recall  that  all 
the  previous  numbers  given— the  sampling  rate,  1000  per  s,  the  quantity 
averaged  25,  and  the  total  time  of  the  sample  10  s— are  all  given  to  identify 
a  region,  not  a  specific  number.  These  numbers  can  all  be  adjusted  to  permit 
512  data  points  to  be  obtained.  Here,  the  two  blocks  of  data  are  available 
with  the  information  of  the  water  vapor  microstructure  necessary  for  deriving 
wind  velocity  value. 

Focusing  the  processing  on  the  Information  of  Interest  in  the  data  is  the  next 
point  to  be  considered.  This  focusing  is  done  by  frequency  filtering  the  two 
averaged  data  blocks.  To  obtain  an  Indication  of  the  frequency  components  of 
interest,  an  example  will  be  used  with  a  10  m/s  wind.  The  highest  frequency 
in  the  data  is  caused  by  the  small  irregularity  passing  through  a  sample 
volume.  As  an  example,  a  4-m  irregularity  moving  along  In  the  wind  will  cause 
2.5  Hz  components  in  the  data.  The  4-m  size  Is  about  the  smallest  change  that 
should  be  considered.  This  comes  about  because  of  the  10  m  aR.  Little  more 
than  1  Irregularity  at  a  time  can  occur  In  the  sample  volume.  The  4-m 
variation  moving  along  in  a  10-m  wind  will  pass  through  the  sampling  point  In 
0.04  s,  hence  2.5  Hz.  As  in  the  next  case,  a  10-m  irregularity  in  the  water 
vapor  microstructure  will  take  1  s  to  pass  through  one  volume,  hence  1  Hz.  As 
the  Irregularity  being  considered  gets  larger,  for  example  up  to  100  m  in 
extent,  the  frequency  component  In  the  data  becomes  0.01  Hz.  With  this  wind 
velocity,  it  would  pass  through  one  volume  In  10  s— the  total  period  of  the 
sample.  One  complete  cycle  of  this  frequenqy  In  the  data  is  of  little 
value.  A  good  point  at  which  to  set  the  low  frequency  cut-off  is  0.02Hz. 
When  these  two  data  blocks  have  been  properly  filtered,  they  are  cross- 
correlated  to  produce  aT.  When  aT  is  obtained,  the  wind  velocity  is  D/aT  *  V. 
V  Is  the  wind  velocity  component  on  a  line  through  the  two  sampling  volumes. 

High  RH  conditions  will  probably  be  a  more  reliable  measurement  than  using 
dust  and  haze  as  the  tracer,*  *•  but  with  this  technique,  information  to 


’Barber  T.  L. ,  and  J.  B.  Mason,  "A  Transit-Time  Lldar  Wind  Measurement,  A 
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derive  a  wind  velocity  by  tracking  dust  and  haze  Irregularities  Is 
available.  Remember  that  using  P^p  and  provides  the  necessary 

Information*  *'  to  derive  the  wind  velocity  from  tracking  natural  dust  and 
haze. 

CONCLUSIONS 

The  feasibility  for  using  the  natural  water  vapor  microstructure  as  a  tracer 
for  remote  wind  measurements  in  a  wet  atmosphere  was  shown.  This 

microstructure  tracking  is  felt  to  be  a  most  appropriate  technique  for 
meteorological  conditions  where  high  RH  is  prevalent. 

Of  interest  here  is  the  relative  variation  in  the  absolute  water  vapor  content 
of  the  atmosphere.  The  variations  or  microstructures  are  measured  at  two 
fixed  points  a  known  distance  apart.  The  wind  velocity  can  be  derived  by 
determining  the  time  of  flight  from  one  measuring  position  to  the  other. 

For  these  procedures  to  produce  the  information  of  interest,  the  DIAL 
measurements  of  the  absolute  water  vapor  content  must  be  taken  with  some  care 
regarding  the  water  vapor  extinction  coefficient  in  aR  sampling  rate,  sampling 
period,  and  data  frequency  filtering. 

When  an  initial  set  of  water  vapor  data  is  collected,  filtering  should  be  done 
carefully  to  focus  the  tracking  process  on  water  vapor  irregularities  in  a 
general  size  range  4  to  100  m.  When  the  processing  is  aimed  toward  tracking  a 
general  size  of  irregularity,  a  reasonable  cross  correlation  is  available; 
hence,  the  wind  velocity  can  be  derived. 

This  technique  has  the  potential  of  working  quite  well  in  certain 
meteorological  conditions.  Recall  that  in  this  same  data,  information  is 
available  to  implement  use  of  water  vapor  or  to  implement  use  of  dust  and  haze 
as  the  tracer.  The  combination  is  certainly  much  more  reliable  in  a  wider 
variety  of  meteorological  conditions  than  either  method  alone. 

The  percentage  of  time  that  adequate  conditions  occur  in  the  water  vapor 
microstructure  to  allow  a  reasonable  wind  measurement  is  an  area  that  should 
be  investigated  further. 
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